Atoms of GROMOS96 Force Field

The derivation of the van der Waals parameters for the aliphatic CH united atoms of the GROMOS96 force field is presented. The parameters n have been adjusted to reproduce the experimental enthalpies of vaporization and vapor pressures or densities of a set of nine alkanes in the liquid state at 298 Ž . K or at the boiling point in the case of methane , using a cutoff radius for the van der Waals interactions of 1.6 nm. Force fields to be used in molecular simulations are bound to the conditions chosen for their parametrization, for example, the temperature, the densities of the systems included in the calibration set, or the cutoff radius used for the nonbonded interactions. Van der Waals parameters for the CH united atoms of earlier GROMOS force fields were n developed using a cutoff radius of 0.8 nm for the van der Waals interactions. Because the van der Waals interaction energy between aliphatic groups separated by distances between 0.8 and 1.4 nm is not negligible at liquid densities, the use of these parameters in combination with longer cutoffs leads to an overestimation of the attractive van der Waals interaction energy. The relevance of this excess attraction depends on the size of the groups that are interacting, as well as on their local densities. Free energies of hydration have been calculated for five alkanes. Q 1998 John Wiley & Sons, Inc. J Comput Chem 19: 535]547, 1998